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1. Introduction 



The circle of ideas which generated Feynman path integrals is contained in works by 
I Dirac and Feynman [U [2] . Especially in the latter work the representation of quantum 
mechanical evolution amplitudes in terms of heuristic integrals on 'path space' was 
developed into an alternative general formulation of quantum dynamics, equivalent to 
the previous formulations by Heisenberg, Schrodinger or Schwinger. 

The miraculous success of Feynman's method in dealing with quantum fields should 
be, however, mitigated by the fact that there is so far no rigorous way to define it in 
terms of conventional measure theory. Namely, the heuristic expression for Feynman 
path integrals is in terms of a complex formal density which does not define a measure. 
Thus the mathematical definition of the objects understood under the name of Feynman 
path integrals posed genuine new problems which have been attacked by different 
methods. For instance, Feynman's path integral in non-relativistic quantum mechanics 
is conventionally viewed as a formal expression which can be given meaning by a specially 
devised limiting process. 

The diverse aspects and approaches make clear that the subject of Feynman path 
integrals should not be considered as a closed one, on the contrary, much work is needed 
on the conceptual mathematical and physical level in order to bring to fruition all 
the beautiful potentialities contained in those ideas. In this direction it was recently 



speculated by G. Svetlichny [31 H] that the secret of the path integral may rest on the 
mutual unbiasedness of temporally proximal bases. He focused the essential problems 
into the following questions: 

(i) For what unitary groups U (t) in L^(R") do the position bases at times and t tend 
to mutual unbiasedness as t — > ? 

(ii) Is there a discrete version of the previous question in a finite-dimensional Hilbert 
space which approximately simulates the propagation of a free particle? 

(iii) What is the information-theoretic nature of the normalization factor A in the short- 
time propagator? 

In this paper we elucidate the true meaning of these questions in the case of a non- 
relativistic quantum particle on the real line. 

In section 2 we recall related notions of complementary observables and of mutually 
unbiased bases. Then the basic concepts of quantum mechanics in finite-dimensional 
Hilbert spaces are introduced in section 3. Section 4 is devoted to the construction 
of finite quantum phase space from the finite Heisenberg group. The group of inner 
automorphisms of finite quantum phase space is described in section 5. After these 
prerequisites, the N x A^-approximation of quantum mechanics on the real line is 
constructed in section 6 and applied to an analogue of a free non-relativistic particle. 
This approximation, being used in the Feynman short-time propagator (section 7) leads 
to the emergence of a Lagrangian as the corresponding local phase and at the same 
time demonstrates mutual unbiasedness of temporally proximal bases (section 8). Our 
derivation also yields the normalization factor A as direct counterpart of the constant 
1 / \/N involved in the definition of mutual unbiasedness. 



2. Complementarity and mutually unbiased bases 

Mutually unbiased bases in Hilbert spaces of finite dimensions are closely related to the 
quantal notion of complementarity. Namely, two observables A and S of a quantum 
system with Hilbert space of finite dimension are called complementary [5], if their 
eigenvalues are non-degenerate and any two normalized eigenvectors \ui) of A and l^j) 
of B satisfy 



\{U^\VJ) \ = 




Then in an eigenstate \ui) of A all eigenvalues bi, . . . ,bN oi B are measured with equal 
probabilities, and vice versa. This means that exact knowledge of the measured value 
of A implies maximal uncertainty to any measured value of B. 

According to W.K. Wootters |6], two orthonormal bases in an A^-dimensional 
complex Hilbert space 

{\u,)\^ = l,2,...,N} and {\vj)\j = 1,2, . . . , N} 

are called mutually unbiased, if inner products between all possible pairs of vectors taken 
from distinct bases have the same magnitude 1 / V^, 

\{u,\vj)\ = for all i,j e {1,2,...,N} . 

V A 

Thus if the system is in the state \ui), then transitions to any of the states \vj) have 
equal probabilities. 
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It is important to note that, in an A^-dimensional Hilbert space, there cannot be 
more than + 1 mutually unbiased bases. It has also been proved that the maximal 
number of + 1 mutually unbiased bases is attained, if is a power of a prime [7] . 

3. Quantum mechanics in finite-dimensional Hilbert spaces 

The mathematical arena for ordinary quantum mechanics is, due to Heisenberg's 
commutation relations, the infinite-dimensional Hilbert space. A useful model for 
quantum mechanics in a Hilbert space of finite dimension A^ is due to H. Weyl 
[H]. Its geometric interpretation as the simplest quantum kinematic on a finite 
discrete configuration space formed by a periodic chain of A^ points was elaborated by 
J. Schwinger [9]. In [10] we proposed a group theoretical formulation of this quantum 
model as well as a finite-dimensional analogue of quantum evolution operator for a free 
particle. 

In an A^-dimensional Hilbert space with orthonormal basis B = {\0),\1), . . .\N — 1)} 
the Weyl pair of unitary operators {Qn, Pn) is defined by the relations 

Qn\p)=ujM, p = 0,l,...,Ar-l, 
Pr,\p) = \p-1 (mod AT)), 

where ujn = exp{2Tri/N) [S] (see also [HI [12] )• If B is the canonical basis of C^, the 
operators Qn and Pn are represented by the matrices 

Qn = diag (l, UN, uJn, ■ ■ ■ , uj^'^) 

and 
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They fulfil a commutation relation 

PnQn = ^nQnPn (1) 

which is analogous to the relation for Weyl's exponential form of Heisenberg's 
commutation relations. Further, = = In, uj^ = 1. 

The finite Heisenberg group is generated hj ujn, Qn and Pn 

= {co'j^Q'nPnI^, J, = 0, 1, 2, . . . , at - 1} . 

It consists of A^^ unitary N x N matrices and is also called the Pauli group. 

The geometrical picture behind the above operators is the following [10]. The cyclic 
group Ziv = {0, 1, . . . A^ — 1} is the configuration space for A^-dimensional quantum 
mechanics. Elements of this periodic chain Zjy provide labels of the vectors of the basis 
B = {\0),\1), . . .\N — 1)} with the physical interpretation that \p) is the (normalized) 
eigenvector of position at p G Z^v. The action of Z^v on Z^v via addition modulo A^ is 
represented by unitary operators U{a) = P^. The action of these discrete translations 
on vectors |p) from basis B is given by 

Uia)\p)=P^\p) = \p-a (mod AT)). 
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The important discrete Fourier transformation is given by the unitary Sylvester 
matrix Sn with elements 

{SN)kp = {p\k) 



/N 

The operator relation 

S^^PnSn = Qn 

shows that the discrete Fourier transform diagonalizes the momentum operator. In other 
words, it performs the transition from the coordinate representation to the momentum 
representation: 

N-l 

\k) = J2\P)iP\k)- 

p=0 



4. Finite quantum phase space 



The following developments will heavily use the finite phase space Tn which is simply 
related to the finite Heisenberg group The center Z{IIn) of the finite Heisenberg 
group is the set of all those elements of IItv which commute with all elements in n^r 

z(nAr) = {(/,o,o)|/ = o,i,...,iv-i}. 

Since the center is a normal subgroup, one can go over to the quotient group Ujy /Z{Iliy). 
Its elements are the cosets labeled by pairs (j, o"), j,a = 0,1, . . . , N — 1. The quotient 
group can be identified with the finite phase space 

Tn = ZnxZn, N = 2,3,.... 

To simplify notation, we denote the cosets corresponding to elements (j, a) of the phase 
space Fat by Q^P'^ without subscripts N 

Q^p'' = {cuj^glp^i / = 0, 1, . . . , iv - 1} . 

The correspondence 

: li^/Zili^) F^ = X : Q^P'^ ^ (j, a), 
is an isomorphism of Abelian groups, since 

= (j, cr) + (/, a') = {j+f,a + a'). 

The group of automorphisms of the quantum phase space F n was studied in jlSl HI] . 
The latter paper considered — instead of cosets — the one-dimensional grading subspaces 
of the Pauli-graded Lie algebra gl{N,C) and studied their transformations under the 
automorphisms of gl{N, C). The subgroup of inner automorphisms was induced by the 
action 

i)x{A) = X~^AX 

of matrices X from GL{N, C). 

Along the same vein we are considering those automorphisms of the above form, 
acting on elements of IItv, which induce permutations of cosets in IIn/Z{IIn)- Operators 
X which induce these automorphisms are unitary. Explicit forms of these operators are 
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given in [13] for prime and in [13] for arbitrary but only for special transformations 
of Tat. 

Automorphisms ip of the given form are equivalent if they define the same 
transformation of cosets in 11^ / Z{U.n)'- 

for all G Z,n x ^Af- The group 11^/2(11^) has two generators, the cosets P and 
Q. Hence if ipy induces a transformation of Il]sf/Z{Il]\[), then there must exist elements 
a,b,c,d E Zjy such that 

y-igr = g«p^ and Y'^PY = Q^P'^. 

It follows that to each equivalence class of automorphisms ipy a quadruple (a, b, c, d) of 
elements in Zjv is assigned. Then we have 

Theorem [14J There is an isomorphism $ between the set of equivalence classes of 
inner automorphisms ipy o,nd the group SL{2,'Z]\f) of 2 x 2 matrices with determinant 
equal to 1 modulo N, 

^{t/jy) = (^^^ a,b,c,deZN; 

the action of these automorphisms on IIn/Z{IIi^) is given by the right action of 
SL{2,Zn) on the phase space = x ^n, 

5. N X N approximation of quantum mechanics on the real line 

An interesting approximation method in quantum mechanics was proposed in E. 
Husstad's PhD. thesis [15j supervised by T. Digernes at NTNU Trondheim. Their 
approach was inspired by an idea of J. Schwinger [16j. 

They approximate quantum operators in for one- dimensional quantum 

systems by A^ x A^ matrices — operators in the Hilbert space /^(Z^r) of finite-dimensional 
quantum mechanics. To this end an auxiliary factor 




is introduced. We have found that for our purpose of approximating the Feynman path 
integral it is still necessary to introduce two additional dimensional quantities: length 
unit a and the corresponding unit of linear momentum h/a. Then the position operator 
is approximated by the multiplication operator in position representation 

QnIp) = ariNp\p). 

For the momentum operator p^ J. Schwinger had the real insight to define it as the 
discrete Fourier transform of q^, implying that p^ has not the form of the generally used 
difference operator. Thus the momentum operator is approximated by the multiplication 
operator in momentum representation 

PN\k) = -r]Nk\k)- 
a 

Schwinger's geometric idea was to identify Zn with a grid in M. For A^ odd, he 
defined a sequence of grids L^v = {ar]]\fp\p = 0, ±1, . . . , ±(A^ — l)/2}. In the limit 
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N ^ oo the grids are becoming denser and at the same time extending to the whole real 
line. The grid serves to embed the finite-dimensional Hilbert space /^(Z^v) isometrically 
in L^(]R, dq) by the map 

_ 1 _i 

where xs denotes the characteristic function of a subset S" C M. The position 
eigenvectors \p) are thus mapped onto narrow normalized wave functions (ppiq) on the 
real line, centered at the grid points and contracting in the limit N ^ oo. Under the 
map X the normalizations of the wave functions \ip) = Sp'^plp) and (Xip){q) are related 
by 

l^pT = ar]N\{1i'){a'nNp)\'^. 



6. Finite-dimensional analogue of quantum free particle 

Finite-dimensional analogue of a quantum free particle was formulated in [10] as a 
discrete Galilean evolution along a finite closed linear chain. The single-step unitary 
time evolution operator Cn proposed there is diagonal in momentum representation 

{j\CN\k) = 5jk^jl'\ 
Transformation to position representation gives 

N-l 



{CN),. = Y.^p\j){j\CN\k){k\a) = ^Y. 



u 

jk ' j=0 

The unitary operator Cn fulfils relations 

Cj^^QnCn = ujnQnPni C^^PnCn = Pn- 

Looking at free evolution in continuous phase space, we arrive to the conclusion 
that operator Cn should be slightly modified. However, first consider the usual one- 
parameter group of unitary operators 

"2 

T(t) = exp(-^^t), tGM, 
h 2m 

describing quantum evolution of a non-relativistic free particle of mass m on the real 
line. The corresponding N x N approximation is 

2 

T^(r) =exp(--^^r£), r G Z, 

where we have introduced a time unit e, since dynamically defined time intervals will 
play a special role. Thus t shall be restricted to integer multiples re, r G Z, of e. The 
time unit e will be chosen so that (in momentum representation) 

TNiT)\j) = ^^P(-^^(^VNjyTe)\j) = u^'^^ r G Z, 

including an additional 1/2 factor in the exponent. [|| Our choice is in agreement with 
a dynamical relation 

ma^ a h 

e = — — or m— = — 

n e a 

\ Non- integer powers oi ujn are understood as complex exponentials tofj = exp(2^u;), w G R. 
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expressing the natural fact that a particle of momentum h/a traverses the distance a in 
time e. Transformation to position representation gives 



TN{r),^ = Y.{p\3){3\TN{r)\k){k\a) = -J^ 



1 -hfr+{p~a)j 



UJ 



N 



jk j=Q 

Now in order to justify the 1/2 factor in the exponent recall that the free time 
evolution in continuous phase space is described by translations along q with constant 
velocity, 

= g(0) + ^t, pit)=piO). 
m 

The corresponding single-step 'translation' in quantum phase space Ftv, 

Q^QP, P 
is equivalent to the following SL{2,Zfyf) transformation 

(1,1) =(1,0) (J ly (0,1) =(0,1) (J I 

Its powers form an Abelian subgroup of SL{2,Zn) isomorphic to Zjv- An easy 
calculation shows that it is implemented by the unitary transformation 

T,,{l)\j)=uj-^^'\j) = Cm\j), 

which, from now on, will be denoted Cm- The modified unitary operator Cjyi now 
fulfils relations 



C^IQnCni — uj^QnPn, C^\PnCni — Pi 



N, 



N 



7. N X N approximation of the Feynman path integral 

Let \q{0), 0) and \q{t),t) be the state vectors of the initial state and of the final state 
of a particle on M at times 0, t, respectively. If S[q] is the classical action functional 
of the particle, the evolution amplitude is according to Feynman formally written as a 
path integral [2] 



{q{t),t\q{0),0) = 




It is understood as a sum over all continuous paths in configuration space. According to 
Feynman's principle of equivalence of trajectories, the contribution of each path should 
have the same absolute value, hence contributes to the sum only a phase factor with 
the phase given by the classical action in units h evaluated along the path. 

In quantum mechanics, the path integral is traditionally defined as a limit via 
discretization based on the division of the time interval, e.g. into n intervals of equal 
duration e = t/n. The evolution amplitude is thus written as a multiple integral [2j 

{qit),t\q{0),0) = 




/ {q{t)\e-^'''\qn_i)dqn-i . . .dqi{qi\e--h'''\q{0)), 
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where qi = q{le) and H is the Hamilton operator. Each factor — the short-time 
propagator — is then identified with an exponential of the short-time action involving 
an approximation of the classical Lagrangian, 

with normalization factor A. For instance, for a non-relativistic particle of mass m 

and one computes (via momentum representation) 

exp(-(p, V{qi))e)dpi = 



2'Kh ^ h e 2m 



— oo 



.2nin£ I i l ,qi+i-qu2 .r, xxx 

'exp(-(-m( ) - V{q^)e), (2) 



I.e. 



Kqi^i^qi) = \<n\— — ^1 -y{qi) and A=(^^^^ 
2 \ e m 



2 



A sequence of g/'s for each ti shall, in the limit, define a path of the system and each of 
the integrals is to be taken over the entire range available to each qi. In other words, 
the multiple integral is taken over all possible paths|§| 

Let us return to our analogue of a free non-relativistic particle. The above approach 
will guide us in our N x N approximation with the short-time propagator induced by 
the unitary operator Cni- In this approximation qi ^ arjNPi, so, for a single time step, 
{qi+i,e\qi, 0) is approximated by 

{qi+ue\qu 0)ariN = {qi+i\e~T^"'^\qi)ar]N = 



{pi+i\Cm\pi) =]^5Z 



1 -i,? 



UJ 



2 



3i+{pi+i-pi)n 



ji=0 



For r time steps 

{pr,Te\po,0) = y ] {Pt\Cni\Pt-i) ■ ■ ■ {Pi\Cni\Po) 



Pi,...,Pt-i 



1 

= {pr\Cm\po) = j;^ 2^^N 

j=0 

The above Gauss-like sum for a single time-step can be summed up using CL. 
Siegel's Reciprocity Formula for generalized Gauss sums [T71 [HI [19] 

|c|-l I |a|-l 

^ ^ni{an'^+bn)/c _ h^^^ni{\ac\~b'^)/ (Aac) ^ ^-■Ki{cn'^+bn)/a 

n=0 " ^ n=0 

§ Recall a quotation from Feynman's thesis (p. 69 of "A point of vagueness is the normalization 
factor, A. No rule has been given to determine it for a given action expression. This question is 
related to the difficult mathematical question as to the conditions under which the limiting process of 
subdividing the time scale, required by equations such as ( 68), actually converges. " 
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valid for a, 6, c G Z, ac 7^ 0, ac + 6 even. Putting a = N with N odd, c = 1, n = ji and 
b = — 2p — 1 with p = pi^i — pi one obtains 



N-l 
31=0 

On the basis of this formula we prefer the operator 

for unitary single-step time evolution. It satisfies simpler relations than Cni, 

Cj^\QnCn2 = QnPni Cj^\PnCn2 = Pn, 
r<-sf)p pj (-<s _ pj+ps 

while inducing the same Abelian subgroup of translations of quantum phase space. With 
operator Cn2 we compute in position representation 

/ \ri \ \ ^ \{pi+i-pi + ^Y 

{Pi+i\Cn2\Pi) = ^=^N 

This result can be interpreted as the emergence of a dimensionless Lagrangian Cn, 

{Pi+i\Cn2\pi) = ^^^"^"'-^""'^ /:^(Pm,P/) = liPi+i-Pi + h'- 

In order to go over to the 1 + 1 space-time and obtain the corresponding local Lagrangian 
Ln, we divide by arj^ and express the short-time propagator 

{qi+i,e\qi, 0) = -^{pi+i\Cn2\Pi) = 
ar]N 



This result can be rewritten 

.2'Kifie 



ViN rn 

where the phase factor appearing in the short-time propagator is seen to define the 
corresponding small increment of the action Ljye which is proportional to the local 
Lagrangian 

_l g^+l-g^ + T ^2 

J^N - ^rn[ ) . 

Note that in the limit ^ 00, we have //at = \J2ti/N and obtain the usual form 
of the short-time propagator for the free quantum particle. 

We close this section with a one-dimensional particle moving in a potential field 
V{q). This potential has been incorporated in the short-time propagator ([2]). To get 
its X approximation, potential V{q) is sampled only at the grid points qi = arj^pi, 
Pi = —{N — l)/2, . . . , (A^ — l)/2. In order to transform V{qi) into dimensionless form it 
should be expressed in the energy unit 



1 h 2 27r h 



used for transforming the kinetic energy to 

V{qi) = V{ariNpi) = -^^^i- 

As a result, potential is represented by a set of N dimensionless constants wi. Thus the 
short-time propagator obtained for a free particle is subject only to a slight modification 
by constants Wi: 

{qi+i,e\qi, 0) = -^{pi+i\Cn2UJn'"'\pi) = 



ar]N ViN 



m \h 2 e 

Also in these formulae the emergence of a local Lagrangian is clearly manifest. 

8. Short-time propagator and mutually unbiased bases 

Let us denote the bases composed of eigenvectors of the operators Q''j^P^ by B^j^a)- 
Unitary operator Cm (or Cni) plays analogous role as operator in our previous 
study of mutually unbiased bases for prime N [20\. There the iterations of generated 
the maximal set of N + 1 mutually unbiased bases 

i^{l,0) — ^ i^(0,l) i^(2,l) — > • • • ^ 1^{N~1,1) 

starting with the canonical basis 'B(i^o)- 

// is prime, identical reasoning as in [20] shows that the iterations of unitary 
operator Cn (indices 1 or 2 omitted) generate in a similar way another maximal set of 
A^ + 1 mutually unbiased bases 

JO '^jv f2 Cjv f2 Cjv f2 Cn Cjv jo 

i^(0,l) ~^ ^{1,0) ^{1,1) ^ ^{1,2) . . . — > 0(i^Ar_i), 

now starting with the momentum basis i3(o,i). Thus the composite unitary operators 
C^S"^^, b = 0, 1, . . . , A^ — 1 produce all the bases of the maximal set when applied to 
the momentum basis. 

In this paper A^ is not restricted to primes but may take arbitrary odd value. 
Notwithstanding this general situation the formulae derived in previous section show 
that the bases appearing in the short-time propagators, i.e. {\p)} = 'S(i,o) and 
{CnU!]^^'-\p)} = B(i^i) are mutually unbiased. Especially, for the short-time propagator 

\{Pi+i\CnujJ;^'\Pi)\ = ^ (3) 



holds and it is this constant absolute value that entails mutual unbiasedness of the 
bases involved. From the physical viewpoint the state evolving after a short time 
interval e carries no information about the preceding state. The trivial information- 
theoretic meaning of mutual unbiasedness therefore consists in the fact that in each 
single evolution step complete loss of information occurs. Physical information is carried 
only by the phase factor whose phase is proportional to the local Lagrangian. 
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9. Conclusions 

The foregoing development of N x N approximations of the Feynman path integral is 
a continuation of our previous studies of quantum mechanics over finite configuration 
spaces [TOl [20] and of discrete path summation [21] . Our approximation method was 
first applied to discrete time evolution of an analogue of a quantum free non-relativistic 
particle in discrete 1 + 1 space-time. The powers of the evolution operator can be 
interpreted as a unitary representation in /^(Z^) of an Abelian subgroup of SL{2,Zfq) 
acting on quantum phase space x Z^. Generalization of the N x N approximation 
of the short-time propagator to a particle moving in a potential turned out to be 
straightforward. 

Summarizing, we have shown explicitly how a Lagrangian arises as a local phase in 
the short-time propagator in the N x N approximation, thus confirming G. Svetlichny's 
conjecture. Our paper also brings definite answers to the questions quoted in the 
Introduction for the case of a one- dimensional non-relativistic particle moving in a 
potential field. 

(i) Our results concern unitary groups U{t) in L^(R) connected with the time evolution 
of a one-dimensional non-relativistic quantum particle moving in a potential field. 
Equation Q shows that position bases involved in the N x N approximation of the 
short-time propagator for arbitrary odd N are mutually unbiased. In this sense also 
in the limit N oo position bases at times and t do tend to mutual unbiasedness 
as t — > 0. 

(ii) Our discrete analogue of quantum evolution of a free particle in a finite-dimensional 
Hilbert space which accurately simulates the Galilean evolution of a free particle 
on the real line was consistently employed throughout the paper. 

(iii) In our paper also the secret of the prefactor in the short-time propagator 
is unveiled. Its dimension must be inverse length in order to compensate the 
integration over q in the Feynman path integral. Hence in the NxN approximation 
it must be related to our unit of length arj^. Due to e = ma'^/h, 



holds for single-step time evolution. This relation shows the trivial information- 
theoretic meaning of A"^: in the dimensionless expression it is a constant 
corresponding to complete loss of information in each single time step. 

In connection with the NxN approximation we would like to point out 
the suggestion of the approximate solution of the continuous Schrodinger equation. 
Namely, Digernes, Husstad and Varadarajan [22] proved a convergence theorem on 
approximation of continuous Weyl systems by x iV Weyl operators Q-'jyP^. Further, 
Digernes, Varadarajan and Varadhan [23] proved a strong theorem on convergence of 
eigenvalues and eigenfunctions of N x N Hamiltonians to solutions of one- dimensional 
Schrodinger equation for potentials satisfying V —>■ +oo as |g| oo, hence possessing 
discrete spectrum. This theorem provides a justification for the approximate solution 
of the continuous Schrodinger equation. Numerical calculations showed that the 
approximation is unexpectedly good even for relatively small values of N. The 
generalization of these results to the case of mixed spectrum remains open. Let us 
note that in quantum optics, discrete phase space Z„ x Z„ is employed in the discrete 
approximation of the quantum phase and the conjugate number operator [21]. 
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